Unified Analytic Electron Interaction Integrals Over Slater Orbitals 

for Diatomic Molecules 



O 

(N 

U 

Or 
< 



E.V.Rothstein 

evrothstein@gmail . com 

A unified formula for the analytical evaluation of two-center ex- 
change, hybrid and coulombic type integrals is presented. 

INTRODUCTION 

A formula for the analytic evaluation of two-center electron interaction 
integrals is presented. The terms in the sum have been checked 1 
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i = orbital number with quantum numbers Ui,li, rrii and screening con- 
stant Si . c — center, centered on nucleus a or b . $j C (j) = Normalized 
Slater Orbital for electron j, in spherical coordinates centered on center 
c. 

I = integral to be evaluated known as the two-center exchange inte- 
gral 

/ = ($ la (l)$36(l)— $2a(2)$46(2)) 
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The S in [Eq. ( 1 )] is the Kronecker delta. The first sum is finite over 
k ( the set of gi k , g 2 K,ri k ,r 2k generated by the four orbitals) . The 
second sum is finite if either or both the /3 's = . If 0i = the sum 
over \x is restricted so the /i + \a\ + g± k is even and the upper limit of 
the summation is g± k + \a\ . If neither of the /3's = the upper limit is 
infinite. Some notation: pj m \cos8) = Associated Legendre function of 
the first kind. ( see Ref. 2 , 8.6.6 ) R = internuclear distance between 
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a and b nuclei. 



®i = —(S 1 + 5 3 ); "2 = — (S2 + S4); fii = — (Si - S3); 

p 2 = ^s 2 -s A y, w = R^ +l mT +l ' 2 

W = 5 1 (a 1 +/3 1 )" 1 - 1 / 2 (a 1 -/3 1 )" s+1 / 2 (« 2 +/3 2 )" 2+1 / 2 (a 2 -/3 2 ) ra4+1 / 2 
W = 4(«i+/3i) ni+1/2 («i-/3i) ri3+1/2 («2+/3 2 ) n2+1/2 (« 2 -/3 2 ) ri4+1/2 
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The square brackets in the upper limit of summation over is \{U — 



or ^(li — | mj| — 1) whichever is an integral. 



9i = Pi + P3 + qi + <?3 + ai + &i + 2ci, 
#2 = + P4 + qi + qA + a-i + b 2 + 2c 2 , 

n = 2(s 1 + s 3 + d 1 ) + q 1 + q 3 -p 1 -p 3 + n 1 + n 3 -l 1 -l 3 -a 1 -b 1 , 

r 2 = 2(s2+S4+d 2 )+q2+q4-P2-P4+n2+n4-l 2 ~h-a2-b2, See footnote 3 
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/3kl+i/2 ' 

= modified spherical Bessel function of the first kind 
( see 10.2.2 Ref. 2 ) 
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If /5 = then, B(ji,g,0, \a\ 
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In practice this latter formula was not used. The ascending 
series formula proved to be more accurate ( less roundoff errors ) 
for j3 7^ and 2k + fj, — \a\ > g, then, 

B{p., g, ft M) = (-J^'^W-. £ 1 ° 

(3) 



-d\ Vl f -d\ r ' 2 f 00 f°° 
A(fi,r 1 ,r 2 ,a 1 ,a 2 , \a\) = [ — J (t^H I I d^d^ 



1 Ji 



x te 2 - i) H/2 (el - i) kl/2 4 a| (6< 2 )glr l (6>i)e- ai?i e-^ 2 



e 



-(a 1 +a 2 ) r/ , I _ I \ i "1 2 2 



22^+1 



(/i+|a|)! 



£h-C-S)(3 



[i£±H] fe fc 



H / W n ^ (fci - rn)! n ^ (fc 2 - n 2 )\ a n 2 

[ l n ( ) + (_)ni+n 2+fcl+fc2 + l e 2(a 1+ a 2)jBi(2ai + ^ 

+ (-) /ci+ni e 2ai E 1 (2a 1 ) + {-) k2+n2 e 2a2 E 1 {2a 2 ) + 7 ] 



&2 1 k2—U2 , , fci 



^ n 2 ^ (fc 2 - ra 2 - j 2 )! 1 (fci - m)! a ? +1 a%> +i * +1 

"2 = 1 J2=0 «i=0 v > I Z 

U I "^Z, 1 on 2 -t-l 

(fci - m)i 1 ] tt(n 2 -t-i)\ < + v 2 2+t+2 

fci! / Ai + n 2 - m - 1\ 1 

( ni ~ ^) ! V fcl " ni / « 2 2+1 («i + a 2 ) fcl+n2 - ni ai 1+1 

^1 (_)fc2+n 2+J2+ l 2 n 2 - t -l / t + kl - Wl \ 1 

+ ^ (n 2 -t-l)! V * jaf 1 (a 1 + a 2 )W+ 1 -«iaf +l)J 
fc 2 ! ^ 1 ^ fcj r 1 



a ni+ii + l a n 2 + l ^ ( ni _ t _l)| V a n+t+2 a n 2+ l 

(") fe2+n2+1 v , 1 e-(^V+H)! ^ fc /> 2/,- 

J J J + 92 M .,| ^ 



X 



X 



(ai + a 2 y+^ +1 a n 2 2+1 2 2 ^! ^ ^ " k 

V 7 K=l \ A* / J= Q \ J / n=0 

2fi-2n \ f/j\ 1 (2/x - 4j - l)(/i - 2j - 1 + |<r|)! 2 2j+1 



- \a\ + k) \nj (2j + l)(n-j)(n-2j-l)\ 

r /x+|<7|-2j-l i 

n /2(//-2j-l-n)\ A-2j-l 



x 

n=0 



£ (-) 



A 4 ~~ 2j — 1 — \a\ J \ n 



fn 1 + f 2 -n 2 
^t^-n.y.^y f 2 -n 2 



"2 



x V ^ 

(n 2 -J 2 )!af +1 K + a 2 )/ 2 +«i-« 2 +i 

A fc 2 ! A fn 2 + fi-rii 

Til , , 



^ (ni - ji)! af + a 2 )/i+« 2 -n 1+ i 

E\{x) = Exponential integral ( 5.1.1. Ref. 2 ) 

7 = Euler's constant ( table 1.1 Ref. 2 ) 

ki — H + \cr\ — 2p + ri, A; 2 = /i + \a\ — 2k + r 2 , 

k{ = n + \a\ - 2k + n, h = fi + \a\ - 2(n + j) - 1 + n, 

/ 2 = A* + kl -2(n + j) - l + r 2 , (4) 
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Hybridlntegral = ($ la (l)$ 3a (l)— $ 2a (2)$ 4b (2)> 
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«i = ft = 2"(5i + ^); «2 = 2"(5 2 + <J 4 ); ft = 2-(<*2 - 

Coulomblntegral = ($ la (l)$ 3a (l)— $ 25 (2)$ 46 (2)> 

n 2 
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«i = ft = — (5i + 5 3 ); «2 = -ft = — (S 2 + S 4 ); 

These changes result in changes in the set of g± k , g 2 k, ri k .r 2k .Ck ■ 



